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Abstract 

We describe the nonzero temperature (T), low frequency (uj) dynamics of the 
order parameter near quantum critical points in two spatial dimensions (d), 
with a special focus on the regime Tiuj <C &bT. For the case of a 'relativistic', 
0(ra)-symmetric, bosonic quantum field theory we show that, for small e = 
3 — d, the dynamics is described by an effective classical model of waves 
with a quartic interaction. We provide analytical and numerical analyses of 
the classical wave model directly in d = 2. We describe the crossover from 
the finite frequency, "amplitude fluctuation", gapped quasiparticle mode in 
the quantum paramagnet (or Mott insulator), to the zero frequency "phase" 
{n > 2) or "domain wall" (n = 1) relaxation mode near the ordered state. 
For static properties, we show how a surprising, duality- like transformation 
allows an exact treatment of the strong-coupling limit for all n. For n = 
2, we compute the universal T dependence of the superfiuid density below 
the Kosterlitz-Thouless temperature, and discuss implications for the high 
temperature superconductors. For n = 3, our computations of the dynamic 
structure factor relate to neutron scattering experiments on Lai.ssSro.isCuO^ 
and to light scattering experiments on double layer quantum Hall systems. 
We expect that closely related effective classical wave models will apply also to 
other quantum critical points in d = 2. Although computations in appendices 
do rely upon technical results on the e-expansion of quantum critical points 
obtained in earlier papers, the physical discussion in the body of the paper is 
self-contained, and can be read without consulting these earlier works. 
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I. INTRODUCTION 



A number of recent experiments have probed the long- wavelength, low frequency, nonzero 
temperature (T) dynamics of the order parameter associated with a T = quantum critical 
point in a two spatial dimensions (d). These experiments include: 

(i) Neutron scattering measurements have mapped out the T, wavevector, and frequency 
dependence of the dynamic spin structure factor in La 2 _ a ;Sr a ,Cu04 for x ~ 0.15 The 
measurements over an order of magnitude in T, and over three orders of magnitude in the 
static susceptibility, are consistent with the presence of a nearby quantum critical point to 
an insulating ordered state with incommensurate spin and charge order ('stripes' 0). 

(ii) Double layer quantum Hall systems at filling factor v = 2 exhibit ground states with 
different types of magnetic order ||. Recent light scattering experiments H have probed 
the fluctuation of the magnetic order parameter in the vicinity of the quantum transitions 
between the states. 

(in) Microwave measurements || of the magnetic penetration depth of high temperature 
superconductors with a number of different T c 's show that the superfluid stiffness satisfies 
the scaling relation 

P ' [T) (i.D 



where ^f p is an apparently universal function. This is precisely the behavior expected in 
the vicinity of a quantum critical point where T c — > ||, such as a superfluid- insulator 
transition. In this paper, we shall provide an explicit computation of the universal scaling 
function ^ p in a model system. We believe our approach and strategy can be generalized to 
models which include some of the additional physics contained in recent discussions P,P|,|TU[| 
of the T dependence of the superfluid density in the high temperature superconductors. 

Motivated by these disparate experimental systems, this paper will present an analysis 
of the long-wavelength, nonzero temperature order-parameter dynamics in the vicinity the 
simplest, interacting quantum critical point in d = 2: that of a 'relativistic', n- component, 
bosonic field <j) a , a — 1 . . . n. However, our ideas and approach are expected to be far more 



general, as we shall discuss further in Section |TV|. The 0(n)-symmetric quantum partition 



function for the field (f) a is given by (in units with h = ks = 1 which we use throughout) 
Z Q = J V<f) a (x,T)exp ^- J d d % J q drCoj 



Cq ~2 



4(^) 2 + (Vx0 a ) 2 + (r c + r)^ 
cr 



+ 1 (0 ■ (") 



Here x is the (i-dimensional spatial co-ordinate, r is imaginary time, c is a velocity, and r c , 
r and u are coupling constants. The co-efficient of the (f) 2 a term (the 'mass' term) has been 
written as r + r c for convenience; we will choose the value of r c so that the quantum critical 
point is precisely at r = 0. So the T = ground state has spontaneous 'magnetic' order for 
r < with (<p a ) ^ 0, and is a quantum paramagnet with complete O(n) symmetry preserved 
for r > 0. The quartic non-linearity proportional to u is relevant about the Gaussian fixed 
point (u — 0) for d < 3, and is responsible for producing a non-trivial quantum critical 
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theory with interacting excitations. Higher-order non-linearities are irrelevant about this 
quantum critical point. 

In addition to being an important and instructive toy model of an interacting quantum 
critical point in d = 2, the field theory ( |1.2|) also has direct applications to experimental 
systems. We will briefly note these now, and discuss them further in Section |TV| . For n — 2, 
Zq describes the transition between superfluid and Mott-insulating states of an interacting 
boson model: + i(p2 is the superfluid order parameter and the quantum paramagnet is a 
Mott insulator. For n = 3, (f) a plays the role of a magnetic order parameter measuring the 
amplitude of the incommensurate, collinear spin density wave in the experiments of Ref. [|TJ . 
The n = 3 case also describes the quantum Hall experiments of Ref. jf], where <ft a now 
measures the difference in the magnetization of the two layers. 

An important tool in the analysis of Zq is the e expansion, where 

e = 3-d (1.3) 

The structure of the e expansion for the T > properties of Zq has already been extensively 
discussed in two previous papers, hereafter referred to as I and II |TT| . We will now 



summarize the main results of these papers, and then turn to a description of the specific 
purpose of this paper. Although the present paper builds upon on these earlier works, an 
attempt has been made to make all of the physical discussion in the body of the paper 
self-contained; earlier technical results are used in the appendices. Some physical results 
from I are summarized in the caption of Fig p], which shall form the basis of our subsequent 
discussion. 

In I, the properties of the phase diagram in Fig[l] were analyzed in an expansion in e. In 
particular, detailed e expansion results were obtained for the dynamic susceptibility 



- f 1/T dr [ d d x f:(0 a (x,r)0 Q (O,O))e-^—), (1.4) 



where k is the wavevector, uj n the imaginary frequency; throughout we will use the symbol 
uj n to refer to imaginary frequencies, while the use of uo will imply the expression has been 
analytically continued to real frequencies. It was found in I that for the static susceptibility 

X (k) = x(k,uj n = 0), (1.5) 

an expansion in powers of \fe held over all regions of the phase diagram of Fig [I], apart from 
a small window in the immediate vicinity of the line of finite temperature phase transitions 
(in this window, the problem reduces to one in classical critical phenomena, and this shall 
not be of interest to us here). So in a sense, the static theory was weakly coupled for 
small e, and this allowed for a satisfactory theoretical treatment of the crossovers in x{k)- 
A completely different situation held for the dynamic properties, and in particular for the 
spectral density lmx{k,u>): the e expansion was found to fail badly, and led to unphysical 
results for small k and u. In particular, in region A, this failure in the computation of the 
dynamic properties appeared for wavevectors smaller than ck ~ yfeT and frequencies smaller 
than uo ~ \ftT . In the e expansion, the low frequency spectral density is given by an integral 
over the phase space for the decay of excitations into multiple excitations at lower energy; 
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however it does not self-consistently include damping in these final states, and this leads 
to unphysical results. In other words, determination of the values of lmx{k,uj), for small 
k and uj requires solution of a strong coupling problem, even for small e, and is dominated 
by the relaxation of excitations with energies of order or smaller than y/eT. Similar results 



hold also for the expansion in 1/n [|12|,|T3j. It is this strong coupling problem which will be 
addressed in this paper. 

Although, as just discussed, the results of I for static properties were adequately com- 
puted at low orders in the e expansion, even they had significant qualitative weaknesses when 
extrapolated to the physically interesting case of e = 1, d — 2, apart from also not being 
quantitatively very accurate. In particular, we know that the line of non-zero temperature 
phase transitions in Fig [l] is not present for n > 3 i.e. T c {r) = for these cases. In contrast, 
leading order e expansion results of I have a T c (r) > for all n. Furthermore, for n = 2, there 
should be a jump in the value of the superfluid density at T c in d = 2, and clearly, this also 
does not appear at any order in the e expansion. We shall address all of these problems in 
this paper, along with the dynamical problem indicated above. We will do this by an exact 
treatment of certain thermal fluctuations directly in d = 2, while the remaining quantum 
and thermal effects (for which the case d = 2 plays no special role) are treated by a low 
order e expansion. We shall claim that this hybrid approach leads to a more quantitatively 
accurate determination of both static and dynamic properties in the high temperature (or 
quantum critical) region A of Fig [l]. Our approach will lead to a computation of the scaling 
function in (|1.1|), containing the universal jump in the superfluid density at T c . 

Before turning to a discussion of our strategy in solving the strong-coupling dynamical 
problem, let us also review the results of II. This paper examined transport of the conserved 
charge associated with the continuous O(n) symmetry of Zq, for n > 2 and e small. The 
T and r dependence of the conductivity was examined using a perturbative expansion in e, 
especially in region A. It was found that for small e, the most important current carrying 
states were bosonic particle excitations with energy Ek ~ T, and momentum k ~ T/c 
(contrast this with the typical energy of order \feT which dominates relaxation of the order 
parameter, as discussed above, and in I). The damping and scattering of the current carrying 
states with Ek ~ T was adequately computed by the e expansion of I, as they were out of 
the region of k, uj space where the weak-coupling expansion broke down. Also, because Ek 
was not much smaller than T, the occupation number of these bosonic modes could not 
be approximated by the classical equipartition value, T/e fc , but required the full function 
l/(e £fe / r — 1) for quantized Bose particles. The transport of charge by these particles was 
analyzed by the solution of a quantum Boltzmann equation in II. All the analysis of II was 
systematic in powers of e, and included only the leading non-trivial terms. The present 
paper will develop a new strong-coupling approach in d = 2, but it will be applied only to 
the low frequency order parameter dynamics; the transport properties of the new approach 
will be examined in a future publication. 



We are now ready to outline the strategy of this paper. We will begin in Section |TA 
by recalling the approach of I for the computation of static properties in e expansion. We 
shall show that a straightforward modification of this approach allows an exact treatment 
of the most singular thermal fluctuations in directly in d = 2, allowing us to obtain results 
which have all the correct qualitative features for all values of n, and are also believed to 
be quantitatively accurate. The low frequency dynamic properties will then be considered 
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in Section |TB 



A. Statics 

The main idea of I was to analyze Zq in two steps. In the first, all modes (j) a (k, u n ) with 
a non-zero frequency, u n ^ were integrated out in a naive e expansion. This produced 
an effective action for the 4> a (k,uj n = 0) modes, which was subsequently analyzed by more 
sophisticated techniques. Our first step here will be identical to that in I; so we define 

,1/T 

$Jx) = T / dT<j>Jx,T). (1.6) 
Jo 

After integrating out the modes with a non-zero u n , the effective action for takes the 

form 



Z = J V3> a (x) exp J d d x C 



2 



1 r 



U /.o\2 



(VAr + fi^+^^j • (1.7) 



The values of the coupling constants R and U will be discussed shortly. We are treating 
the consequences of the non-zero uo n modes at the one loop level, and at this order the co- 
efficient of the spatial gradient term does not get renormalized. This approximation means 
that effects associated with wavefunction renormalization and the quantum critical exponent 
rj have been neglected: this is quite reasonable, as n takes rather small values at the 2 + 1 
dimensional quantum critical point. Also, these two loop effects were considered at length 
in I, and were found to be quite unimportant. 

We note in passing that Z in (|T77) is designed to apply as a model of Zq in (|1.2|) in region 



A of Fig |T|. It also contains the initial crossovers into regions B and C, but there are some 
subtleties as T — > in regions B and C. As we shall see in Appendix B, the limits T — > 
and e — > do not commute: the leading physics for very small T can be properly captured, 
but there are some subleading effects which are not accounted for in an approach based on 
the e expansion. These caveats also apply to the dynamics to be discussed in Section 1TB. 



For our remaining discussion, it is crucial to understand the properties of Z as a contin- 
uum, classical field theory in its own right. Our strategy here will be obtain the universal 
properties of this continuum theory directly in d = 2. Actually (|1.7|) does not define the 
theory Z completely, as some short distance regularization is needed to remove the ultravi- 
olet divergences. A priori, it might seem that there is no arbitrariness in choosing the short 
distance regularization, as it is uniquely provided by the underlying quantum theory Zq. 
However, as will become clear now, it is actually possible to choose a 'virtual' short distance 
regularization of Z at our convenience, provided we properly match certain renormalized 
couplings with those obtained from the true quantum regularization due to Zq] we will work 
with a 'virtual' lattice regularization of Z here. So, what are the short distance singularities 
of Z ? From standard field theoretic computations ||14|| , it is known that for d < 3, the 



model Z has only one ultraviolet divergence coming from the 'tadpole' graph shown in Fig|2| 
(there are some additional divergences, associated with composite operators, which appear 
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when two or more field operators approach each other in space: we will not be concerned 
with these here — however, these are important in a consideration of transport properties, 
and will be discussed in a future paper). So all short scale dependence can be removed 
simply by defining a renormalized coupling R related to the bare coupling R by 

~ „ rT1TT fn + 2\ f 1 / 11 d d k 1 . , 

R = R-TU[——) / - — - — . 1.8 

V 6 J J (27i) d k 2 + R y 1 

Here the expression l/(k 2 +R) should be read as a schematic for the low momentum behavior 
of the propagator. At higher momenta of order 1/a (for lattice regularization, we will 
choose a to be the lattice spacing), the propagator can have a rather different momentum 
dependence and this has to be accounted for in computing the integral in (|1.8|). Also notice 
that we have performed the subtraction with a propagator carrying the renormalized 'mass' 
R. For d > 2, this is not crucial and we can equally well define the subtraction with 
a massless propagator 1/k 2 : this procedure was followed in I, and has the advantage of 
leading to a simple linear relation between R and R. Here, we are interested in working 
directly in d = 2, and then such a massless subtraction would lead to an infrared divergence. 
So we are forced to perform the subtraction as in ( |1.8| ). Indeed in d = 2, ( |1.8| ) evaluates to 

R = R-TU — ^- )ln(C/Ra 2 ), d = 2 (1.9) 

where C is a regularization dependent, non- universal constant. This clearly shows that it is 
not possible to set R = in the subtraction term. We also note an important property of 
( p. ,8| , |I~9|) which is special to d — 2. Clearly, we have assumed above that R > 0. However 
( |1.9[ ) shows that the bare mass R ranges from — oo to +oo as R increases from to oo. So 
it is no restriction to consider only positive values of R, as that allows us to scan the bare 
mass in Z over all possible negative and positive values (this is not true for d > 2 as the 
reader can easily check from ( p..8| ): then we do need values of R < 0, while defining the 
renormalization with a massless propagator, to access all values of R). In particular, as we 
will show, we will be able to access both the magnetically ordered and disordered phases 
of Z for R > in d = 2. We can also interpret Q1.9Q in a renormalization group sense as 
defining the scale-dependent effective mass R at a length scale a, in a theory with a fixed, 
positive R: so even in a theory with R > 0, it is possible to have a significant window of 



length scale, where the scale-dependent mass R is less than zero. We will see in Sections [D I 
and |Iy| that this interpretation is very helpful in understanding the origin of 'pseudo-gap' 
physics in the quantum critical region. 

After short distance dependencies have been removed by the simple renormalization in 
( |1.9|) , all correlators of Z are expected to be finite in the limit a — > 0, and are universal 
functions of the renormalized couplings R and U. Actually, instead of working with R and 
U, we shall find it more convenient to use R, and the dimensionless Ginzburg parameter, Q, 
defined by 

TV 

G = (L10) 

as our two independent couplings; the ratio Q gives an estimate of the strength of the 
non-linear fluctuations about the mean-field plus Gaussian fluctuation treatment of Z. So, 
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provided we express everything in terms of R and Q (and not R and U), the properties of 
Z are regularization independent and universal functions of R and Q. 

At this point, the correct approach towards computing the static properties of the un- 
derlying quantum problem Zq should be quite clear. 

(i) First, we compute the values of the effective couplings R and Q, defined in ( |1.7|) , ( |1.8|) 
and (|1.10| ), by integrating out the nonzero uo n modes in flL2|) . This is carried out by an 



extension of the approach developed in I, and our results are presented in Appendix B. In 
the most interesting high T region A of Fig [l], these couplings take the following values to 
leading order in e 



R = e 



n + 2\ 2n 2 (T/c) 



n + 8J 3 



= *-, 4W5 ■ (LID 

y/2(n + 2)(n + 8) 

As one moves out of the region A, these couplings become smooth, monotonic, and universal 
functions of r /T l l zu . In particular, Q obeys the simple scaling form 



S = *4^i7^J, (1.12) 

where ^lg is a universal function, and C\ is a non-universal scale which can be eliminated if 
the argument of is related to the actual T = energy gap or spin stiffness (as discussed 
in I and Appendix B). As one moves in Fig [1] from Region C to A to B, the argument 
of tyg decreases uniformly from +oo to — oo, while the value of Q increases monotonically 
and analytically. The value of <&g(0) is given in fll.ll ). Similar results hold for R, which 



decreases monotonically from region C to A to B. In the body of this paper, we shall regard 
R and Q as known functions of r and T which can be looked up in Appendix B and I. 
(u) Then, we study Z directly in d = 2. We do this mostly with our own convenient choice 
of a (virtual) lattice regularization, but we will be careful to express all physical response 
functions in terms of R and Q after using ( |1.8|) ; we will explicitly show that our results 
become independent of a, for small a, when this is done. Finally, for these couplings R and 
Q, we use the Zq imposed values given in Appendix B and (|1.11| ), to determine the true 
physical response functions. The values of R and Q vary non-trivially as a function of r and 
T, and the final results then contain the crossovers between the different regions of Fig |l| 

Simple, engineering dimensional analysis shows that the static susceptibility of Z, and 
therefore also of Zq, has the form (as shown in I) 

T X (k) = (*«(*)*«(-*)> = ^ S\ , (1-13) 

where ^(k,Q) is a universal function of its two arguments. One of our primary tasks here 
shall be the computation of this universal function directly in d = 2. For small Q, this can 
be done in naive perturbation theory; as is clear from ( |1.11| ) such a perturbation theory is 
applicable in region A for small e, as Q ~ \ft. However, in d = 2 and e = 1, the value of 
Q in fll.lip is not particularly small, and the resulting perturbation theory for the static 



7 



susceptibility is not quantitatively reliable. We shall present the results of straightforward 
numerical simulations carried out on a small workstation, which give a reasonably accurate 
determination of the universal function \1/ in d = 2, except when Q is extremely large. 
Somewhat remarkably, precisely in d = 2, we shall also be able to make exact statements in 
the strong coupling limit Q — > oo through a duality-like transformation, and this will provide 
a useful supplement to the numerical results. In summary, by a combination of weak-coupling 
perturbation theory, an exact strong-coupling 'duality' mapping, and numerical simulations, 
we shall obtain fairly complete knowledge of ^ directly in d = 2. 

In d — 2, for the cases n = 1,2, there are critical values Q = Q c where Z exhibits 
phase transitions (in the universality classes of d = 2 classical Ising and Kosterlitz-Thouless 
respectively), which appear as singularities of the function \1> at Q = Q c . The values of Q c 
will be determined numerically (for n = 1, Q c was obtained in Ref [|1^]). For Zq, these phase 



transitions, of course, reflect those along the T > full line within region B in Fig |l|. 

An important property characterizing the low temperature phase for d = 2, n = 2 and 
Q > Q c is the spin stiffness, p s {T). Simple dimensional arguments similar to those leading 



to ( |1.13[ ) now show that its temperature dependence obeys 

p s {T) = T^ p {Q) (1.14) 

where ^ p is a universal scaling function; it is closely related to the experimentally measurable 
\l/ p in ( |1 . lp , and this will be discussed in Section [TVI. Clearly, 4f p vanishes for Q < Q c . 



B. Dynamics 

As we have emphasized in I and II, the dynamical properties of region A of Fig |l] in d = 2 
are especially interesting because they are characterized by a phase coherence time, and an 
inelastic scattering time, which are both universal numbers times T. Consequently, thermal 
and quantum fluctuations play equal roles in the dynamical theory; this novel regime of 
dynamics was dubbed quantum relaxational in Ref [jnj. Here we argue that for the case 
where e is small, and for the long wavelength relaxational dynamics of the order parameter 
only, it is possible to disentangle the quantum and classical thermal effects. The central 
reason for this is that the predominant modes contributing to the relaxation of the order 
parameter fluctuations at long wavelengths have an energy of order 

c\fR~^T, (1.15) 



from ( |1.11|) . It must be emphasized that these modes carry negligible amounts of current, 



and the transport properties continue to be dominated by excitations with energy of order T 
even for small e, as we have discussed in II. As the energy in (|1.15| ) is parametrically smaller 



than T, the occupation number of the typical order-parameter modes is 

1 T 1 , 

(1.16) 



The second term above is the classical equipartition value. So we can conclude that there is 
an effective classical non-linear wave model which describes the long wavelength relaxation of 
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the (p a fluctuations. Quantum effects then appear only in determining the coupling constants 
of this effective classical dynamics. 

So what is the classical wave model describing the relaxation of the order parameter 
modes ? To leading order in e, the required model can be deduced immediately from some 
simple general arguments. First, we have the important constraint that the equal time 



correlations must be identical to those implied by Z in ( |1.7|) . Second, to endow the <3> Q field 
with an equation of motion, we clearly need to introduce a conjugate momentum variable U a . 
The kinetic energy associated with this momentum is clearly given by the time derivative 
term in Lq in ( p..2[ ). Furthermore, we know that the co-efficient of this gradient term is not 
renormalized at order e when the high frequency modes are integrated out. So we assert 
that the required dynamical model is specified by the following partition function over a 
classical phase space 



Z c = J V® a (x)VTl a (x)exp 



n 

T 



n 



J d d x [l [c'Ul + (V,$ Q ) 2 +R$l} + | ($9 2 } . (1.17) 



Notice that we can perform the Gaussian integral over ITq, exactly, and we are left then 
with the original co-ordinate space in ( |1.7| ): this is the usual situation in classical statistical 
mechanics, where momenta decouple from the static analysis. To complete the specification 
of the classical dynamical model, we need to supplement (|1.17| ) with equations of motion. 
These are obtained simply be replacing the quantum commutators associated with classical 
Poisson brackets. So we have 

{$ a {x), X\ p {x')} pB = 5 aP 5{x - x'). (1.18) 

The deterministic, real time equations of motion are then the Hamilton- Jacobi equations of 
the Hamiltonian 7i, and they are given by 

<9<I> 

= c 2 n a , (i.i9) 

and 



dt 



{n Q (x),H} 



P.B. 



v 2 x $ a -m a -^($l)$ a , (i.2o) 



The equations ( |1 . 1 7|) , ( |1.19D and Q1.20D define the central dynamical non-linear wave model 



of interest in this paper. We are interested in the correlations of the field $ a at unequal 
times, averaged over the set of initial conditions specified by ( |1 . 1 7|) . Notice all the thermal 
'noise' arises only in the random set of initial conditions. The subsequent time evolution 
obeys Hamiltonian dynamics, is completely deterministic, and precisely conserves energy, 
momentum, and total O(n) charge. This should be contrasted with the classical dynam- 
ical models studied in the theory of dynamic critical phenomena [|l(|[17]], where there are 



statistical noise terms and an explicit damping coefficient in the equations of motion. 
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The dynamical model above has been defined in the continuum, and so we need to 
consider the nature of its short distance singularities. Our primary assertion is that for 
d < 3, the only short distance singularities are those already present in the equal time 
correlations analyzed in Section f^]. These were removed by the simple renormalization in 



( |1.8| ), and we maintain this is also sufficient to define the continuum limit of the unequal time 
correlations. These assertions rely on our experience with the structure of the perturbation 
theory in Q presented in I, and on the consistency of the numerical data we shall present 
with the scaling structure we describe below. 

Assuming that introducing R as in ( ]1.8[ ) allows to take the limit a — > 0, we can deduce 
the scaling form of unequal time correlations by simple dimensional arguments. We define 
the dynamic structure factor, S(k,u) by 

1 r°° r n 

S(k,u) = - dt d d x ^($ Q (x,t)$ ct (0,0))e- i{fc ^ ) . (1.21) 
n J-oo J a=1 

Notice that, unlike ( |1.4| ), this involves an integral over real time, t. Comparing with the 
equal time correlator in (|1.13|) , we clearly have the relation 

/OO rj[.) 
^W) s{k,u) ' (L22) 

However, what is the relationship between S(k, uo) and the physically appropriate quantum 
dynamic susceptibility x{k, obtained by analytically continuing ( |1.4| ) ? By analogy with 
( |1.21| ) we can define the physical, quantum dynamic structure factor 

1 r°° r " 

S Q (k,u) = - dt d d x ^(0 a (x,t)0 a (O,O))e-^-^, (1.23) 
n J-oo J a=1 

where it is understood the time evolution is now due to the quantum Hamiltonian implied 
by Zq, and so is the thermal average. This structure factor obeys an exact fluctuation 
dissipation relation to x(k,u) defined in flO]): 

S Q( k ^) = 1 _ 2 e - u /T hR x(k,u). (1.24) 

We can relate the dynamic structure factors Sq and S only under the conditions that the 
dominant spectral weight of excitations is at an energy smaller than T. As argued earlier, this 
is the case here by (|1.15|) for e small. Assuming this condition we have S(k,uj) w Sq(k,u) 
and 

2T 

S(k,u) w — Jmx(k,u). (1.25) 

Finally, we can write down the scaling form obeyed by S(k,u); from the arguments of 
the previous paragraph, and simple engineering dimensional considerations as in ( |1.13| ), we 
obtain 

= « „), ( , 26) 
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where *f>sc{k, oo, Q) is a dimensional universal function is an even function of 10. The prefactor 
of ([OH) has been chosen to ensure that this function has a constant integral over frequency 

00 duo — 

— y Sc (k,co,g) = l, (1.27) 

-00 Z7T 

as follows immediately from (|1.22p , or from ( |1.25| ) after use of the Kramers- Kronig represen- 
tation of x(k) in terms of lmx{k, 00). We will obtain information on the structure of ty$c for 



d = 2 in Section |TJ. We shall be especially in interested in the to — > limit of S(k, to) which 
describes the long time correlations of the order parameter $ Q ; this limit is not accessible 
in perturbation theory even for small Q, as was shown in I. 

The outline of the remainder of the paper is as follows. We will present details of our 
analysis of the static properties of Z in Section |I[ The main achievements of this section are 
of a technical and quantitative nature, and there are no qualitatively new physical results; 
some readers may wish to skip this section and go directly to Section |TTT| . In Section |TJ we 
will present our numerical results on the long time dynamics of the model Zc- A synthesis 
of our results in the context of its experimental implications will then appear in Section [IV]. 
Some technical details, including a summary of needed results from I and II are presented 
in the Appendices. 

II. STATICS IN TWO DIMENSIONS 



This section will examine the classical model (|1.7| ) directly in d — 2. We will obtain 
essentially complete information on the static susceptibility, x(/c), by a combination of weak 
coupling (Section [II A| ), strong coupling (Section [II B|) , and numerical methods (Section [II C| ). 



The exact duality-like transformation will be described in Section II B 



A. Weak coupling 

For small Q, we can perform a familiar Feynman graph expansion in the quartic coupling 
in Z. At order Q 2 we obtain at k = 



" 2 ' N j(i, i, i)0 a + o(0 8 ; 



(2.1) 



where J(l, 1, 1) is a number defined in ( |C19| ), and ( j2.1|) is clearly consistent with the scaling 
form ( p. .13) . The order Q 2 correction starts becoming important for Q « 15, which then is 
roughly where the crossover to strong coupling occurs. 

The spatial correlations of $ decay exponentially on a scale of order 1/y/R. Neglecting 
the order Q 2 corrections, we get from a Fourier transform of the l/(k 2 + R) propagator at 
large \x\ 

1 n T 

- £<<^)<*>a(0)> = e-N/g, (2.2) 

n a=i y8n\x\/^ 

where the correlation length £ is 
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1 + 0{G 2 ) . (2.3) 



B. Strong coupling 

We will now consider correlators of Z in the limit Q —>■ oo. Quite remarkably, exact 
information can be obtained in this limit too. The key is an ingenious proposal made some 
time ago by Chang ]TS|] for the n = 1 Ising case, but which appears to have been forgotten 
since. Chang proposed a strong-to-weak coupling mapping for n = 1, which has the flavor 
to a duality transformation. Here we will review his mapping, and show that closely related 
methods can be applied to all n. 



The argument begins by noting that in the limit R — > 0, the bare mass R in ( |1.8|) tends 
to — oo. So at short scales, the effective potential controlling the fluctuations of \Q a \ will 
have a negative curvature at the origin and a minimum at a non-zero value of \& a \. This 
suggests that we renormalize the theory with a negative renormalized co-efficient of the 
term. So we replace (|1.8|) by 



n- Rd TIj( n + 2 \ * k 1 (2 A) 

R --—- TU {—)J (2n)*k* + R D > (2 ' 4) 

where we have introduced a new renormalized 'dual' mass Rd > 0; the factor of 1/2 in the 
first term on the right hand side is for future convenience, and has no particular significance. 
While both (|2.4|) and ( |1.8|) have a non-universal cutoff-dependence in their momentum in- 
tegral, this disappears when we combine them to eliminate R and obtain 

2 V 24vr J R K J 

This equation can be solved to yield Rd > as a function of U and R. As we will discuss 
shortly, such a solution exists only for R small enough. 

Assuming the existence of a solution fl2.5|) for Rd, we now have a new renormalized 
theory in which the local effective potential for <3> Q fluctuations has the form 

-^ + T(*S) a + -"> ( 2 - 6 ) 



where the ellipses represent counter-terms arising from ( |2.4| ) which will cancel the cutoff 
dependencies order by order in U. So if this renormalized theory is weakly coupled, will 



fluctuate around a non-zero mean value of order ^3Rd/U. But is this the case ? Clearly, 
this depends upon the value of a 'dual' dimensionless coupling, Q D , which is the analog of 

TU 

Qd = -j^- (2.7) 

What is the value of Qd ? This can be obtained by combining (|2.5| ), (|1.10|) and ( |2.7| ) into 
the dimensionless equation 
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which can be solved to yield as a function of Q. A plot of Qd versus Q is shown in Fig |3|. 
A straightforward analysis of ( |2.8| ) shows that there is no real solution for Qd for Q < 54.2756 
for n — 1, for Q < 40.7069 for n = 2, and Q < 32.5655 for n = 3. For Q larger than these 
values, there is a solution for Qd which decreases monotonically from (^o) max = 127r/(n + 2) 
as Q increases, and as Q — > oo it obeys 

So Qd is small as Q becomes very large, and this 'dual' problem is therefore in a weak 
coupling limit, as we had hoped. 

We have so far been expressing all universal physical properties of ( |1.7|) in terms of Q 
and R. However, for large Q, it is clear that we can freely trade these couplings for Qd and 
Rd'- we define Qd from (|2.8|) and relate Rd and R by the relationship 



RdQd = RQ, (2.10) 

obtained by comparing (|1 . 10|) and (|2.7|) . We will mainly use Rd and as our independent 
couplings in the remainder of this subsection. 

It now remains to do a weak coupling analysis in powers of Qd- This is straightforward 
for n=l, but the cases with continuous symmetry, n > 2, have to be treated with some 
care. All of these analyses have been discussed in Appendix C, and we will present the final 
results for different value of n in the following subsections. 



1. n = 1 



For n — 1, we have from ( |C5|) for large \x 
($ a (x)$ Q (0)) = iV 2 + 

where the correlation length 



T 



1 



R 



D 



and the spontaneous magnetization was computed in (|C8| ): 



iV n 




1 - 



56? 



D 
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J(1,1,1) + 0(Q D ) 



(2.11) 



(2.12) 



(2.13) 



The numerical constant J(l, 1, 1) also appeared in (|2.1|) , and its value is given in ( |C19|) . So 
the correlator approaches A^q exponentially on the scale £. In momentum space, the static 
susceptibility has the form 



X ( i ) = M (ar)VW + _^ i 



12.14) 



which is of the form ( |1.13|) . The presence of the delta function indicates true long-range 
order which breaks the $ — > — $ symmetry for small Qd (large Q). 
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2. n = 2 



For n = 2, from (|C9|) and the arguments below it, we can conclude that for |x| 3> 1/ V-^d, 
there is a power-law decay in the order parameter correlator 



t (* a {x)*a{0)) = (l + ^^Od + 0{Gl) 



a=l 



\x\\ Rj 



(2.15) 



where 7 is Euler's constant, and the continuously varying exponent, 77, is related to p s (T), 
the exact renormalized spin stiffness towards 0(2) rotations, by 



T 



' 2np s (T)- 

We computed p s (T) in Appendix |C 2| in a perturbation theory in Q D and found 

p s (T) 3 Q D 



T 



ID 



36 



(2.16) 



(2-17) 



which is consistent with the scaling form ( 1.14 ). The power-law decay in ( |2.15| ) corresponds 
to a quasi-long-range XY order in the two-component planar order parameter $ 2 ). In 
momentum space, the quasi-long range order implies a power-law singularity in the static 
susceptibility at k = 0: 



X(k) 



3r(l-7j/2) 



7T 



2^(77/2) g D Rfh*- 



(2.18) 



which is consistent with the scaling form (|1.13|) . 



3. n > 3 



For n > 3, no long-range or quasi-long-range order is possible. Correlations always decay 
exponentially at sufficiently long scales, but the correlation length does become very large 
in the strong coupling limit. We conclude from the analysis in Section |C 3| that the ultimate 
long-distance decay of the correlation function has the form 



n T p 

E(^w^(o))=c 1 ^r i)/{ "- 2)e 

VD 



where C\ is a universal number, and the correlation length is given by 

l/(n-2) / 



(2-19) 



1 -r / " 



D 



n-2 



e(n - 2)Q D 
48^ 



exp 



6tt 



{n-2)Q D 



(2.20) 



The static susceptibility can be deduced from the results of Ref [ff9|,|13j, and is given by 



X(k) = C 2 



^xV(n-2) 



(2.21) 
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for small k, where C2 is a universal number (C2 ~ 1.06 for iV = 3 ||20||), / is a smooth scaling 
function considered in Refs fll9lJl3l with /(0) = 1. Notice that, unlike the cases n — 1,2, 
there is no singularity in x{k) at k = for small Gd (large Q); instead x(0) becomes 
exponentially large (~ £ 2 ), has an exponentially small width (~ l/£) in momentum space, 
but remains a smooth function of k. 



C. Numerical results 

We now have an understanding of the properties of Z both in the limits Q — > and 
G — > 00 in d = 2. For small G, we have the result ( p.2| , [2.3| ) showing the <3> Q correlations 
decay exponentially in space due to the fluctuations of n $ Q modes about $ a = 0, while 
fl2.1| ) shows that the static susceptibility x(0) ~ 1/R. For large Q we have the results for 
the static susceptibility in (|2.14j ), Q2.18| ), and (|2.21| ). We will examine the manner in which 
the system interpolates between these limits in the following subsections. 



1. n = 1 



The delta function in (|2.14|) indicates the presence of long-range order for sufficiently 
large Q. This delta function is expected first appears at a critical value G = Gc by a phase 
transition in the universality class of the d = 2 classical Ising model. The value of Q c was 
determined numerically in a recent Monte Carlo simulation in Ref []15||, which found 



61.44. 



(2.22) 



2. n 



The quasi-long range order implied by ( [2.15D and (|2.18|) is expected to be present for 



G > Gc, and to vanish at Gc by a Kosterlitz-Thouless transition. The low temperature phase 
is characterized by the spin stiffness which obeys the scaling form in ( |1.14| ). The behavior 



of the scaling function \I/ p for large G (small Gd) was obtained earlier in ( p. 17] ). We expect 
that 

^ P (G < Gc) = 0, (2.23) 
while precisely at Gc its takes the value specified by the Nelson-Kosterlitz jump 

V P (G = Gc) = ~ (2.24) 

7T 

We performed Monte Carlo simulations to obtain more information on the functional 
form of and the value of Gc- We discretized (|1.7|) on a square lattice of spacing a, and 
used an L x L lattice with periodic boundary conditions. The Monte Carlo sweeps consisted 
of two alternating steps. First we updated both the amplitude and phase of ($1 + 1^2) on 
each site by a heat bath algorithm |^TJ. Then we applied the Wolff cluster algorithm [22 
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to rotate the phase of sites on clusters by a random angle. As we are interested in fairly 
large values of Q, where *Sf p is non-zero in the thermodynamic limit, it is more appropriate 
to use the dual couplings, Qd and Rd in testing for the appearance of the continuum limit. 
In particular, we need Rdo 2 <C 1, and we used values around R^a 2 ~ 0.04; this was found 
to yield a-independent susceptibilities, as we will display more explicitly in our discussion 
of the n = 3 case. 

Our numerical results for p s {T) are shown in Fig [|. The stiffness was measured by 
evaluating the expectation value of the appropriate current-current correlation function im- 
plied by the Kubo formula. The results are presented by plotting (p s (T) /T) j '(3 / 'Q d) versus 
and we will now discuss the reason for this choice. From ( |2.17| ) and our discussion in 
Appendix B we see that Qd vanishes linearly as T — > 0, and that 

lim ^ = -4v (2-25) 

where we have now emphasized that ^ is a function of T, as was also noted in the scaling 
form ( |1.12|) ; this relationship guarantees that vertical co-ordinate in Fig || becomes unity as 



Qd — ► 0. Further, if we approximate the T dependence of Qd by Qd(T) ~ 3T/p s (0) (this 
relationship is not exactly true), then the vertical axis in Fig |] becomes p s (T)/p s (0), while 
the horizontal axis is T/p s (0). 

Notice that for small Qd, the results for p s (T) are approximately independent of L, while 
they become strongly L dependent around Qd ~ 3, as would be expected in the vicinity of 
the Nelson-Kosterlitz jump, which is present only in the infinite L limit. We can make quite 



a precise estimate of the position of this jump by fitting [23fl the L dependence of p s to the 



following theoretically predicted |24]] finite-size scaling form 



m = U( 1 + —1—) ; (2.26) 

A and Lq are free parameters, determined by optimizing the fit. The best fit values of 
A are shown in Fig |5|. The value of Qd at which A = 1 determines the position of the 
Kosterlitz-Thouless transition, and in this manner we determine 

Q Dc = 2.747. (2.27) 

Finally, using Q2.8Q we get 

Q c = 102. (2.28) 



3. n = 3 

There is now no phase transition as a function Q, and the susceptibility exhibits a smooth 
crossover from the weak-coupling form (|2.1|) to the strong-coupling limit (|2.21| ). We obtained 
numerical results for x(0) at intermediate values of Q, and the results are shown in Fig|B|. A 
range of values of L and Ra 2 were used, and the excellent collapse of these measurements in 
Fig |^ indicates that we are studying £ in ( |1.7| ) in the continuum and infinite volume limits. 
The weak-coupling prediction of ( |2.1| ) is also shown, and this is seen to work only for very 
small values of Q. 
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III. DYNAMICS IN TWO DIMENSIONS 



Finally, we turn to the central problem of dynamic correlations. We generated initial con- 
ditions for $ a as described in Section [II C[ for Il a by the simple independent Gaussian distri- 
butions specified by ( |1.17D , and then integrated the equations of motion (|1.19|) and (|1.20| ) by 
a fourth-order predictor-corrector algorithm. Correlations of $ Q at unequal times were then 
measured, and in this manner we obtained the correlation function / d 2 x (<& a (x, t)<& a (Q, 0)). 
The results are shown in Figs [7|-|5] for n = 1, 2 and 3 respectively, for a series of values of 
Q. The values of Q were chosen to be around the quantum-critical value evaluated 
directly in e = 1. Also in Fig |9|, we show results for different values of Ra 2 , and their inde- 
pendence on this parameter is evidence that we are measuring the universal values in the 
continuum limit. 

There is a simple, and important, trend in the dynamical correlations with increasing Q. 
For small Q, the k = correlations show a clear damped oscillation in time. These oscillations 
represent amplitude fluctuations in $ Q about a minimum in the effective potential at $ Q = 0. 
The damping of the oscillations increases with increasing Q, until the oscillations disappear 
entirely for Q large enough. 

The Fourier transform of the data in Figs to frequency directly gives us the dynamic 
structure factor, and the scaling function $c defined in ( |1.26| ). The results for this are shown 
in Figs |10^[P2^ for n = 1, 2 and 3 respectively. Notice that S(0, 0) is clearly always non-zero. 
Consequently, by the fluctuation-dissipation theorem, fll.24| ) or ( |1.25[ ), lmx(0,u;) ~ ui for 
small to. The perturbative computations in I did not obey this simple and important low 
frequency limit, and so this sickness has been cured by the present non-perturbative, but 
numerical, computation. 

The small Q regime of amplitude fluctuations discussed above in the time domain, trans- 
lates now into a peak in S(0, u) at finite frequency of order ~ c\J~R. As Q is reduced, we 
move out of the high T region A in Fig [l], and into low T region C on the quantum paramag- 
netic side. This finite frequency, amplitude fluctuation peak connects smoothly with a sharp 
peak associated with a quasiparticle excitation of the quantum paramagnet. Of course, once 
we are region C, the amplitude and width of the peak can no longer be computed by the 
present quasi-classical wave description, and we need an approach which treats the excited 
particles quasi-classically. 

Now consider the opposite trend of increasing Q towards the low T region B on the 
magnetically ordered side. As Q increases, the peak broadens and eventually, the maximum 
moves down to zero frequency. For n > 2, it is natural to interpret this dominance of low 
frequency relaxation as due to "phase" or "angular" fluctuations of $ a along the contour of 
zero energy deformations at a fixed non-zero | $ a | . Of course the fully renormalized effective 
potential necessary has a minimum at $q, = because this is a region without long range 
order; nevertheless, there must be a significant intermediate length scale over which the 
local effective potential has a minimum at a non-zero value of and the predominant 
fluctuations of $ Q are angular. This can also be seen from the relation ( |1.9| ): crudely, we 
can imagine varying a at fixed R to determine the effective mass R on a length scale a — we 
see that for large Q, there is a significant scale over which R is negative, which prefers a 
locally non-zero value of \<f) a \, and allows for low-energy phase fluctuations. 

We can understand the nature of the dynamics in the limit Q — ► oo by arguments 
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analogous to those made for the statics. We will restrict our attention here to n = 3, and 
other cases are similar. We saw in Section II B 3 that the statics where described by the 
d = 2 0(3) non-linear cr-model. In a similar manner we can argue that the dynamics will 
be given by the dynamical extension of this model considered by Tyc et al. ||20|| , and the 
three-argument universal scaling function \l/ $ c in ( |1.26| ) will collapse to the two-argument 
universal scaling functions of Tyc et al. in the limit Q — > oo. This collapse is similar to 
the transformation of (1.13) to ( |2.21p in the same limit for the statics. Consistent with the 
interpretation in the previous paragraph, the description of the Q — > oo limit of the dynamics 
given by the model of Tyc et al is described by a model in which \<p a \ is constrained to have 
a fixed length. Further the results of Tyc et al show a large uj = peak in S(0,uj) |20| . 
which is consistent with the trends observed here with increasing Q. 

The above description of the origin of the low frequency relaxation in the continuum high 
T region A due to angular fluctuations in <p a clearly relies on the existence of a continuous 
symmetry for n > 2. However, closely related arguments can also be made for n = 1 by 
appealing to the low-energy mode arising from the motion of domain walls between ordered 
regions with opposite orientations. 

We have implicitly assumed above that 'amplitude' and 'angular' fluctuations are mu- 
tually exclusive phenomena, but this is clearly not true in principle. Even in a region with 
angular fluctuations, there can be an amplitude mode involving fluctuations in \<f> a \ about 
its local potential minimum. Such a situation would be manifested by a simultaneous peak 
in 5(0, u) both at uj = and at a finite frequency. It is apparent from Figs |10l - |T2| that such a 
situation never arises in a clear-cut manner. Apparently, once angular fluctuations appear, 
the non-linear couplings between the modes is strong enough in d = 2 to reduce the spectral 
weight in the amplitude mode to a small amount. However the amplitude mode does not 
completely disappear — there is a clearly visible shoulder in the n = 1 Fig [TO] for Q = 35, 
indicating concomitant angular and amplitude fluctuations. These results on the difficulty 
of observing an amplitude mode for large Q in d = 2 (low T in region B) connect smoothly 
with the T = response of the magnetically ordered state of the quantum theory Zq — the 
latter is reviewed in Appendix A, and we find there that the amplitude mode is swallowed 
up in the spin-wave continuum for the continuous symmetry case. 

For the quantum-critical region A in Fig [l], we should use the value of Q in ( l.llj) . At 



e = 1, this result evaluates to Q = 35.5 for n = 2, Q = 29.2 for n — 2, and to Q = 24.9 
for n = 3. If we take this value of Q seriously, then we see from Figs |IT)|-|r2| that all 
cases are quite close to the border between amplitude and phase fluctuations, when the 
peak in S(0, uj) moves from non-zero to zero frequency. Amplitude fluctuations are however 
somewhat stronger for n = 3 (when there is a well-defined peak at a non-zero frequency), 
while angular/do main- wall relaxational dynamics is stronger for n — 1 (when there is a 
prominent peak at u — 0). 

There is a passing resemblance between the above crossover in dynamical properties as a 



function of Q, and a well-studied phenomenon in dissipative quantum mechanics [25,26,27 



the crossover from 'coherent oscillation' to 'incoherent relaxation' in a two-level system cou- 
pled to a heat bath . However, here we do not rely on an arbitrary heat bath of linear 
oscillators, and the relaxational dynamics emerges on its own from the underlying Hamil- 
tonian dynamics of an interacting many-body, quantum system. Our description of the 
crossover has been carried out in the context of a quasi-classical model, but, as we noted 
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earlier, the 'coherent' peak connects smoothly to the quasiparticle peak in region C of Fig|l|; 
in this latter region the wave oscillations get quantized into discrete lumps which must then 
be described by a 'dual' quasi-classical particle picture. 



IV. IMPLICATIONS FOR EXPERIMENTS 



We first summarize the main theoretical results of this paper. It is convenient to do this 
in two steps: first for the statics, and then the dynamics. 

For static properties, we presented a rather complete analysis of the classical field theory, 
Z in ( |1.7| ), of a n component scalar field $ a with a quartic self interaction. We motivated the 
study of Z here as an effective theory of the static fluctuations of the quantum model Zq, in 
( |1 . 2| ) , but it is clear that Z has a much wider domain of applicability. The static properties 
of almost any quantum model in d = 2 with an 0(n)-symmetric order parameter should be 
satisfactorily modeled by Z. Of course, the T-dependent values of the coupling constants, R 
and Q will then be different, and depend upon the specific underlying quantum model. So, 
for instance, if one of the phases was a c?-wave superconductor, and had gapless fermionic 
excitations, then the subleading corrections to expressions for the couplings R and Q in 
( |B19[ ) would change. Apart from a traditional weak coupling analysis of Z, we introduced a 
surprising, exact solution of the strong coupling limit by a duality-like transformation. We 
also interpolated between these limits by Monte Carlo simulations. Among our main new 
results was a computation of the T dependence of the spin stiffness, p s (T) for n = 2, and 
shown in Fig f|. When combined a knowledge of the temperature dependence of Q as in 
the scaling form ( |1.12| ) it leads to a prediction for p s (T) consistent with the form ( |1.1| ). As 
a first pass, we can combine the approximate low T prediction which follows from ( |B19| ), 
Qd ~ 3T/p s (0), with Fig|] and obtain an explicit prediction for the function ^ p in ( |1 . 1| ) . 

Our studies of dynamic properties were somewhat more specialized to the quantum model 
Zq, although we expect that closely related methods can be applied to other models. Our 
approach relied heavily on the specific value of the 'mass' R obtained in the high T limit 
of Zq — we used the fact that c\fRjT ~ y/e < 1 to argue for an existence of a dynamical 
model of classical waves. This model is defined by the ensemble of initial conditions (|1.17|) 
and the equations of motion ( |1.19|) and (|1.20|) . The universal dynamical properties of this 
model were studied by numerical simulations directly in d — 2, and results are summarized 
in Figs [FHT2|. As we pass from region C to A to B in Fig [l], the value of the coupling Q 
increases monotonically. The dynamics shows a crossover from a finite frequency, "amplitude 
fluctuation", gapped quasiparticle mode to a zero frequency "phase" {n > 2) or "domain 
wall" (n = 1) relaxation mode during this increase in Q. 

One of the primary applications of our results is to the superfluid-insulator transition 
for the case n = 2. Our approach offers a precise and well-defined method to describe 
the physics of strongly fluctuating superfluids above their critical temperature, with an 
appreciable density of vortices present. Near a quantum critical point, our results show that 
such superfluids (region A of Fig [I]) have a reasonably well-defined order parameter, with 
non-zero over a significant intermediate length scale. Strong phase fluctuations [§] are 
eventually responsible for the disappearance of true long-range order. The evidence for this 
picture comes from our dynamic simulations, which show a well-formed relaxational peak 
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at zero frequency in the dynamic structure factor. The trends in the dynamic structure 
factor as a function of Q then support the interpretation that this peak arises from phase 
relaxation. We speculate that our results can be extended to deduce consequences for the 
electron photo-emission spectrum of the high temperature superconductors, along the lines 
of Refs P|L0||. We imagine, in a Born-Oppenheimer picture, that the fermionic quasiparticles 
are moving in a quasi-static background of the $ Q field. Then the photoemission cross section 
can be related to a suitable convolution of the electron Green's function and the dynamic 
structure factor of Under such circumstances, we believe that a zero frequency, phase 
relaxation peak in the dynamic structure factor of $ Q will translate into a weak 'pseudo-gap' 
in the fermion spectrum. 

A separate application of our results is to zero temperature magnetic disordering tran- 
sitions for the case n = 3. Recent neutron scattering measurements of Aeppli et al. |T|] on 
La 2 _ x Sr x CuO"4 at x = 0.15 are consistent with quantum-critical scaling with dynamic critical 
exponent z = 1 and anomalous field exponent rj ~ 0, suggesting proximity to an insulating 
state with incommensurate, collinear spin and charge ordering (the collinear spin-ordering 
ensures that a single n = 3 vector order parameter is adequate; coplanar ordering would 
require a more complex order parameter and is not expected to have rj close to 0). Their 
measurements have so far mainly focused on the momentum dependence of the structure 
factor, and are well fit by a Lorentzian squared form. We have not computed such momen- 
tum dependence in our simulations here, but general experience with scaling functions in 



d = 1 suggests that such a form is to be expected in the high T regime A [p8[ . It would be 



quite interesting to examine the u dependence of the structure factor in future experiments, 
and compare them with our results in Fig [12|. As we discussed in Section |TT|, it is the smaller 
values of Q, which have a non-zero frequency peak in S(0,u) in Fig [12|, and which lead to 
a 'pseudo-gap' in the spin excitation spectrum. Compare this with our discussion of the 
n = 2 superfluid-insulator transition above, where a zero frequency peak in 5(0, oS) at larger 
values of Q was argued to lead to a fermion pseudo-gap. It is satisfying to note that if we 
are to observe both a spin and a fermion 'pseudo-gap' in the experiments, then the trend in 
the required values of Q with n is consistent with fll.ll ). 



Continuing our discussion of the application of n = 3 to the high temperature super- 
conductors, it is also appears worthwhile to remind the reader of nature of the magnetic 
spectrum in the magnetically disordered side, in region C of Fig [1] [|29|JT3|1 . Here there is 
a sharp triplet particle excitation above the spin gap, which is only weakly damped. It is 
interesting that such excitations have been observed at low T in YBa 2 Cu06+ x |30[| . How- 
ever, their eventual interpretation must await more detailed studies and comparison with 
the situation in La 2 _ a ;Sr :r Cu04. 

It is tempting to combine our discussion of n = 2 and n = 3 models above for the 
high temperature superconductors, into a single n = 5 model |31| . However, there does not 
appear to be any reason for the resulting theory to be even approximately 0(5) invariant. 
Moreover, we have seen above that we need the freedom to independently vary Q for the 
n = 2 and n = 3 subsystems, at moderately large values, to obtain a proper description of 
the physics; for n = 5, the value of Q in (|1.11|) is quite small, and would lead to a rather 
sharp gap-like structure in the dynamic structure factor of the n = 5 order parameter in the 
high T limit. 

Finally, we mention recent light scattering experiments H on double layer quantum Hall 
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systems which have explored both sides of a magnetic ordering transition. Simulations at 
n = 3, but in the presence of a magnetic field can lead to specific predictions for this system. 
The experimental results appear to show the appearance of a spin pseudo-gap at high T, 
which is consistent with our results so for n = 3. 
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APPENDIX A: SPECTRUM OF THE ORDERED STATE OF THE QUANTUM 

THEORY AT T = 



We address here the issue of amplitude fluctuations of \<p a \ in a state with magnetic long- 
range order. We will do this by examining the response functions of the quantum theory Zq 
in ( IO| ) at T = 0. These were computed in I by the e expansion — for small e, there is a well- 
defined peak in the spectral density of the longitudinal response functions, corresponding 
to amplitude oscillations of |0 a | about a non-zero value. Indeed, these results were used 



by Normand and Rice |32| to argue that such an amplitude mode will be observable in the 
insulator LaCu02.5, which is believed to be near a d = 3 quantum critical point. 

We will consider the case d = 2 here, using the large n expansion. Unlike d = 3, we will 
find here that such an amplitude mode is not visible in d = 2 because the cross-section for 
decay into multiple spin-wave excitations is too large. This result also connects smoothly 
with our T > 0, d = 2 simulations in Section pTI , where again we found little sign of such 
an amplitude mode. 

Large n results for the two-point <p a correlator in the direction parallel (x\\{k,u)) and or- 
thogonal (x±{k, u)) to the spontaneous magnetization were given in Appendix D of Ref |]T3" 



They are universal functions of p s (0), c, and the ground state spontaneous magnetization, 
N , and in particular, they do not depend upon whether the underlying degrees of freedom 
are soft spins (as in Zq) or vectors of unit length (as in Ref To leading order in 1/n, 

the results are 



N 2 



X\\(k, 



p.(0) [k 2 - (u/c) 2 ] 
N 2 



Ps(0) ^k 2 - (co/c) 2 [^k 2 - (uj/c) 2 + 16p s (0)/cn 



(Al) 



So 



as expected, the transverse correlator has a simple pole at the spin-wave frequency, 
oj = ck. On the other hand, the longitudinal correlator only has a branch-cut at u = ck — the 
spectral density vanishes for u < ck, and decreases monotonically for u > ck. In particular, 
there is no pole-like structure at a frequency of order p s (0)/n, which is the expected position 
of the amplitude mode. 
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APPENDIX B: COMPUTATION OF R AND U 



This appendix discuss the values of R and U obtained by integrating out the non-zero 
imaginary frequency modes from Zq. Such a computations was already discussed in I, but 
here we will present the modifications necessary due to the slightly different renormalization 
used in flLSp . We will also present new computations within the magnetically ordered state 
in region B of Fig [l]. 

First, in Section [B 1| , we will follow the paramagnetic approach of I, which computes 
parameters for r > 0, and then extrapolates to r < by a method of analytic continuation 
in r, which is valid for T > 0. This method works without any hitches in regions A and 
C of Fig |T[ In principle, it is also expected to be valid within all of region B, but the 
results becomes progressively poorer as the limits T — > and e — > do not commute for 
r < 0. We will then present, Section [B 2| an alternate computation which begins within the 
magnetically ordered state of region B and then directly computes the 'dual' couplings Rd 
and U . 



1. Paramagnetic approach 



We begin by noting that the value of the bare coupling r c in ( P-.2| ) is 

2\ r d d k r duo 1 



-u 



n 



6 



(2n) d J (2tt) k 2 + (u/c) 



(Bl) 



to leading order in u. We assume that r > 0, and so it is valid to integrate out fluctuations 
in cf) a about <p a = 0. Integrating out the u n ^ modes from Zq in this manner to leading 
order in u, and comparing resulting effective action with Z, we find straightforwardly 



R = r + u 
U = u — u 



n + 2 



2 i n + 



6 



d d k 

(27l) d 

d d k 



1 



duo 



T So k2 + Mc) 2 + r J (2ir) k 2 + (co/c) 2 
1 



(2n) dT ^ (k 2 + (co n /c) 2 + r) 2 



(B2) 



The remaining task is, in principle, straightforward: we have to combine (|B2j ) with 
and evaluate the resulting expressions to obtain our final results for R and U. In the vicinity 
of the quantum critical point in Fig [I], the resulting expressions should be universal functions 
only of T, c, and an energy scale measuring the deviation of the ground state couplings from 
the r = point; for this energy scale we choose either the ground state spin stiffness, p s (0) 
(for r < and n > 1), or the ground state energy gap (which is A_ for r < 0, n — 1, and is 
A + for r > and all n). 

One reasonable approach at this point is to solve (|B2|) and (|1.8|) for R and U directly in 
d = 2 by a numerical method. This will give results for R and Q which are valid everywhere 
in the phase diagram of Fig |l| — the resulting value of Q will increase monotonically from 
to oo as one moves from region C to A to B. Moreover, R will remain positive everywhere. 
However, the results will not be explicitly universal and will depend upon microscopic pa- 
rameters from the theory Zq — the values of u and the momentum cutoff A. 
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Universality of the final result can only be established order by order in e. In the 
remainder of this subsection we will evaluate ( |B2|) in such an expansion in e. This surely 
reduces the accuracy of our final estimates for R and U ; it should not be forgotten, however, 
that we subsequently study Z directly in d = 2, and so the scaling functions in (|1.13|) and 
( |1.26| ), whose arguments are related to R and U, are known much more accurately. We will 
find that our leading order in e result for R vanishes when T becomes sufficiently small in 
region B — then the e expansion can no longer be considered adequate for estimating R in 
d = 2. The vanishing of R is acceptable for small e (d > 2) however, because we can then do 
a massless subtraction in ( (L.8|) and negative values of R merely place the system well within 
the magnetically ordered state. An alternative, e-expansion computation of R for systems 
well within region B will appear in the following subsection. 

The techniques for reducing (fB^) and ( |1.8| ) into a universal form in the e expansion have 
been discussed at length in I, and also in Ref @. Using these methods we find to leading 
order in e 



c 2 R + 2vre ( ) cTVR 



cU 



6e(T/c)( 3 - d > 
(n + S)S d+1 



c 2 r 



1 + e 



1 + e 



(20 + 2n - n 2 
2(n + 8) 2 




+ eT 2 



n 



n + 8 



G 



c 2 r 



(B3) 



where S d = 2/[r(d/2)(4 7 r) a! / 2 ] a phase space factor, v = 1/2 + e(n + 2)/(4(n + 8)) to this 
order in the e expansion, and fi is a short-distance momentum scale which can be eliminated 
by re-expressing r in terms of physical energy scales. The function G(y) was given in (D8) 
of II: 



G{y) 



-2 / dq 



In 2q 



2 (cosh(yg 2 + y) 



y 



q 2 + y 



2 vV + 1/e. 



(B4) 



This form of G(y) is valid for both negative and positive y (when the argument of the square 
root is negative we use the identity coshzx = cosx) and is easily shown to be analytic at 
y = where 



G(0) 



2tt 5 



dG/dy(0) = 2.453808582 .... 



(B5) 



The simple form G(y) ~ G(0) + y(dG/dy(0)) is actually a reasonable approximation for 
G(y) over a wide range of values of y. For y — > oo, we can show from (|B4j ) that 



G(y 



oo 



ylny 



+ 2-K^fy- + J8n^e 



-y/v. 



(B6) 



The expression (|B3|) for U is identical to that in I, while that for R differs only in that 
the term proportional to T\/R~ on the left hand side was absent in I. This difference is of 
course due to the subtraction term with mass R in ( p..8|) . Because of the presence of this 
term, ( |B"3| ) is reasonable only as long as there is a solution with R > 0. For low enough T 
in region B there will be no such solution, and then the present method breaks down as a 
method for estimating R in d = 2; an alternative approach will be discussed in Section [B~2 
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We complete this appendix by relating r to physical energy scales to leading order in e. 
The reader can easily verify that when we eliminate r in (|B3|) by the following expressions, 
the arbitrary scale /i disappears to the appropriate order in e. The following relations can 



also be used in to similarly eliminate r from the expressions in Section [B 2 



For r > 0, the ground state has an energy gap, A + for all n. Then, from I, we have 

(B7) 





For r < 0, we have to distinguish n — 1 and n > 2. For n — 1, there is an energy gap, 
A_, and 

(B8) 

For n > 2 it is convenient to use the parameter p s obtained from the ground state spin 
stiffness, which has the engineering dimensions of energy in all d (in d — 2 it is simply 
proportional to p s (0)): 

d = c (d-2)/(d-i) ( 2e Ps(0) Y /(d 1] (B ) 

Ps \(n + 8)S d+1 ) ■ 

Then 

/ 12 + n-2n 2 \ ( -2rY ^ . 

Ps = cp[l-e Al , oN9 Hr • ( B1 °) 



2. Magnetically ordered approach 



Here we will directly compute the 'dual' couplings Rd and U by matching the effective 
potential for $ Q fluctuations in Z in (|1.7|) with that obtained by integrating out the nonzero 
Matsubara frequency modes in Zq in (|1.2|) . We will initially assume that we are in the 



magnetically ordered state, and so r < and R < 0. The effective potential of ( |1.7|) has a 
minimum at |$ Q | 2 = —QR/U, and its curvature at this minimum is —2R/T. We compute 
precisely the same quantities from the free energy obtained from Zq after integrating out 
the nonzero Matsubara frequency modes at the one loop level, and so obtain expressions for 



R/U and R/T. We solve these for R and U, and obtain results which replace (|B2 



R = r + u 



d d k 
(2vr) d 



n + 2 



k 2 + (uo n /c) 2 - 5r (n - l)(k 2 + (w n /c) 2 - r) 



-u 



2(k 2 + K/c) 2 - 2r) 2 
d d k r du 1 



+ 



6(& 2 + (w n /c) 



2\2 



6 



17 



u — u 



d d k 
(2tt) c 



(27r) d 7 2vr A; 2 + (cu/c) 2 
3 

+ 



n 



(/c 2 + (uj n /c) 2 - 2r) 2 6(k 2 + K/c) 2 ) 2 



(Bll) 
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We now use (|2.4j) and evaluate these expressions by the same methods which led to ( |B3| ). 
This replaces ( p3|) by 



c 2 R 



D 



2-r ( H±|') cTjR D = -c 2 r 



cU 



where 



2 \n + 8 

6e(T/c)( 3 " d ) 



(n + 8)S d+1 



1 + e 



(20 + 2n - n 2 
2(n + 8) 2 




= 3G(2y) + (n - 1)G(0) - 9yG\2y) 
G 2 (y) = 9G'(2y) + (n-l)G'(0). 



in 



i)yG'(o) 



(B13) 



There expressions can now be combined with ( |B~8| ) and ( |B10| ) to obtain universal expressions 
for Rjj and U within region B, and also into the crossover into region A. The resulting values, 

53f) because 



when combined with (|2.5|) will not agree precisely with those obtained from 
the computations in the present appendix are good to leading order in e, while the relation 
( [2.5|) is only valid in d — 2. 

Let us look at the values of Rn and U obtained from ( [B12| ) in the limit T — > with 



r < 0. After using (|B6|) it is straightforward to show that 



R 



D 



-2r 



1 + 



9e 



+ e 



n 



2{n + 8) \n + 
ireT /An + 5 
I n + 8 



cU 



6e(-2r) 



c 

(3-d)/2 




V^27- 



[n 



e 



(20 + 2n-n 2 ) 
2(n + 8) 2 
2c 2 r 



2(n + 8) 



[n 



Din 



2^2 



+ 9 + 2(n - 1)G'(0) + 



18ttT 



+ 



(B14) 



where only exponentially small terms of order e~ c ^~^ //T have been omitted. The ln(T) terms 
above are dangerous as they do not have a finite limit as T — > 0. For n = 1, a glance at ( |B14| ) 
shows that all such terms do indeed cancel; so the result ([B12|) remains valid everywhere in 
region B, including in the limit T — > 0. Further it can be checked (using results in I) that 
Qd(T -> 0) = 3T/iV 2 , which agrees with (|T3|). 

However, for n > 2, there is no cancellation of such terms. This indicates a breakdown 
of the e expansion for these values of n as T — ► 0. The low- lying excitations for these 
cases are gapless spin waves. In the present approach, by matching the quantum theory to 
the classical action Z, we are effectively treating these spin waves as classical up to a high 
energy cutoff of cy/—2r, which is the mass of the amplitude mode. In reality, however, such 
spin waves are only classical up to an energy k^T and this is responsible for the 

breakdown of ([B12|) . 

In the remainder of this subsection, we abandon the e expansion, and use some reasonable 
physical criteria (suggested by the above discussion) to estimate Qd and R as T — > 0. Such 
estimates are, by nature, unsystematic, and there does not appear to be any clear-cut 
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procedure by which they can be improved to extend to higher T. We will match the exact 
large Q form of the correlators of Z in Sections [II B 2 and II B 3 with known exact results 
for the low T correlators of the quantum theory Zq. These quantum correlators can be 
universally expressed in terms of p s (0), c, and the ground state spontaneous magnetization 
No 

For N = 2, the predominant fluctuations at low T are angular fluctuations about some 
locally ordered state. So we write 



$i + i$ 2 ~ e . 

The quantum action controlling fluctuations of 9 is 

Ps{0) 



(v x ey + (d T 



(B15) 



(B16) 



Evaluating the two-point correlator from (|B15 ) and ( |B16| ), and fixing the missing normal- 
ization in QB15|) by matching to the actual value of N , we obtain 



£($a(aO$a(0)) =iV 2 exp 



a=l 



d 2 k 



p s (0)J (2tt 



,:ikx 



- 1 



+ 



du 



k 2 + (tu n /c) 2 J 2n k 2 + (u/c) 2 



(B17) 



The expression (|B17| ) is free of both ultra-violet and infra-red divergences, and we obtain 
for large |x| 



^($ Q (x)$ Q (0))=iV 2 exp 



o=l 



T 



27rp s (0) 



(ln{T\x\/c) - In 2 + 7) 



(B18) 



where 7 is Euler's constant. Matching this correlator with that of the classical theory Z in 
(p9|), we obtain two important, and exact results 

w 

£ D (T^0) = — + 



Ps(0) 
T 2 



+ 



(B19) 



where the ellipses represent unknown terms which are higher order in T. 

A similar computation can be carried out for n > 3. In this case, comparing the classical 
result ( |(J29| ) with the quantum results of Refs |T9| , |I3"|| , we find precisely the relations ( |B19[ ). 
So ( gig) holds for all n > 2. 



APPENDIX C: STATIC CORRELATIONS IN THE 'DUAL' THEORY 

This Appendix will show how we can compute correlators of the theory ( |1.7| ) in the limit 
that the Ginzburg parameter Q becomes very large. This is done in the 'dual' formulation 
discussed in Section |11B| : is a theory characterized by a Ginzburg parameter which 
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becomes small when Q becomes large. All of the results will be expressed in a perturbation 
theory in Q D . 

We first perform a naive computation of the O (^-invariant two-point $ Q correlator in 
the original static theory (|1.7| ), assuming R < 0. We will then renormalize it using Q2.4|) , 
and will find, as expected, that the result is finite in the limit a — > 0. The interpretation of 
the result will be straightforward in the for the case n = 1 and will lead directly to ( [2.11|) . 
We will then turn to further computations for the cases n = 2, and n > 3 in subsequent 
subsections. An important ingredient in the interpretation of the results for these cases shall 
be the computation of the change in the free energy of the theory ( |1.7| ) in the presence of an 
external field which couples to the generator of 0(n) rotations: this will allow us to compute 
the renormalized spin stiffness for n = 2, and the precise correlation length for n > 3. 

At the mean-field level, for R < (which is assumed throughout this appendix), $ a 



fluctuates around the average value \$> a \ = \/6\R\/U. So we write 



/61-Rl ~ 

= y-^ a ,i + $<*(>), (CI) 

where $ Q represents the fluctuations about a mean-field magnetization in the a = 1 direc- 
tion. Inserting ( |C1| ) in ( |1 . 7| ) we find that the action changes to 



1 n r 9 i \R\Tl n TT ( n \ 

C =nH (V A) + <W2|ip 2 J + V^P^l £ ^ + 57 £ &a) (C2) 



2 



This action will form the basis for a perturbation theory in [/ (or equivalently by ( |2.7| ), in 
Qd) in the rest of this appendix. Notice that there is a cubic term in this action, and so 

7^ 0. Indeed, the expression for the incipient 'spontaneous magnetization' iVo = 
correct to zeroth order in U, can be obtain in a straightforward perturbative calculation 



from (|C2j) , and we obtain: 



N 



'6\R\ 



TU /■!/« d 2 k ( 1 n 

l--7~7/ 7^77 — -7-7 + 



4|i2|y (2vr) 2 Vjfc2 + 2| J R| 3fc 2 



(C3) 



The fluctuation corrections above come from directions longitudinal and transverse to the 
incipient magnetization respectively. Notice that the latter lead to an infrared divergence for 
n > 1, which correctly suggest that the transverse fluctuations are so strong that the actual 
spontaneous magnetization is 0. For now we will ignore this infrared divergence, as we will 
see that this divergence disappears upon insertion of ( |C3j ) into combinations which are 0(n) 
invariant. To complete the computation of the 0(n)-invariant two point $ a correlator, we 
need the correlations of $ Q . It is sufficient to compute these at tree level, and we have the 
formal expression 

~ ~ v r 1 /® H 2 h p ikx 

<S>Jx)® a (0))=T , (C4) 

1 ' K J (2vr) 2 p + 2|/2|5 ail ' 1 ' 

where there is no implied sum over a on the left hand side. We can now combine (|C1J ), 
(P3|), and (|U4|) to obtain an expression for J2a=i(^oi(x)^ a (0)) correct to zeroth order in U. 
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In this expression we express R in terms of Rd using Q2.4|) , and also expand the resulting 
expression consistently to zeroth order in U. Finally, using the definition oi Qd in (|2.7| ), we 
obtain one of the main results of this Appendix: 



£<*«(*)*«(0)> 
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3T 
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1 + 



Qd f d 2 k 



Akx 



(2 



(n-1) 



e ikx _|_ n 



1 



7T 



A: 2 



k 2 + R D 



0{Q_ 



(C5) 



It is simple to check the crucial property that (05) is free of both ultraviolet and infrared 
divergences; the former arises because the renormalization in ( |2.4|) is expected to cure all 
short distance problems, and the latter because the left hand side is 0(n)-invariant. 

The remaining analysis is quite different for different values of n, and we will consider 
the cases n — 1, n — 2 and n > 3 in the following subsections. 



1. n = 1 



The expression ( |C3|) for the spontaneous magnetization is free of infra-red divergences 
for this case; it can also be made ultra-violet finite by expressing R in terms of Rd using 
(PI). This procedure leads to the simple result that 



iV n = k\ — 



1 + 0{Q%) 



(C6) 



The expression (C5) shows that the two-point $ Q correlator approaches A^q exponentially 
over a length scale 1 / y/Ro- 

It is not difficult to extend (|C6|) to obtain the two-loop 0(Q D ) correction by evaluating 
in a perturbation theory in U under the action (|C2] ). We verified that by expressing 
R in terms of Rd using ( |2.4| ), and expanding consistently to the needed order in U, all the 
ultraviolet divergent terms in the resulting expression cancelled, and we obtained 



Nn 
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The integrals in ( |U7|) are easily evaluated to yield 
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/ 3T 
V yo 
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where the numerical constant J(l, 1, 1) is given later in (C19). 



(C7) 



2. n 



Evaluating 



for x ^> 1/ \[R~d and n = 2 we obtain 
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n <r K r T 

E(^(x)$ Q (o)) = — 



Qd 
6tt 



\u(\x\JR d ) -In2 + T + 



(C9) 



where 7 is Euler's constant. We are now in the magnetically ordered state of the n = 2 XY 
model, and we know this system has power-law decay of spatial correlations. Therefore we 
assert that ( |C9"D exponentiates into the result ( |2.15|) , and identify the spin stiffness p s (T) as 

Mp- = -2- + o(jS o ). (cio) 

1 yD 



We will now compute the next two terms in the small Qd expansion for p s (T) in 



Rather than obtaining these by computing higher order corrections to the correlator in (|C9p , 
we will compute the stiffness directly by examining the response of the free energy density 
to an external field H which couples to the generator of 0(2) rotations. This is equivalent 
to computing the change in the free energy in the presence of twisted boundary conditions. 
In particular, we modify the action in (|L7|) by the substitution 

(V,$ Q ) 2 (V^i - iH§ 2 ) 2 + (V,$ 2 + iH®x)\ (Cll) 

and compute the small H dependence of the free energy In 2(H). The stiffness is defined 
by 

1 Z(H)_ Ps (T) 

v ln Y{o)-^r H + ••• (C12) 

where V is the volume of the system, and the ellipses represent terms higher order in H. 
We compute this expression by first modifying (|Cl| ) to account for the H dependence of the 
mean-field magnetization: 



and expanding the resulting action in powers of H. In this manner it is not hard to show 
that to order H 2 



VH 2 Z(0) TU \ 1V ;/ \ 24|i?| \ U ^ V ' 
-i J d 2 y((® x {y)V y ®2{y) - ®2{y)V y ®i{y)) ■ ($i(x)V,$ 2 (x) - ® 2 {x)V ^{xj)) , (C14) 



where all expectation values are to be evaluated under the action (|C2j) . It is a straight- 
forward, but lengthy, exercise to compute the right hand side of ( p!4j) in a perturbation 
theory in U, which requires enumerating all Feynman graphs to two loops. After formal 
expressions for the graphs have been obtained, we perform the substitution Q2.4p to replace 
R by Ro, and again collect terms to order U 1 . We will not present the details of this here, 
but will state the result obtained without any further manipulations on the expressions for 
the individual Feynman graphs: 
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(C15) 



It is directly apparent that all terms in ( U15|) are individually ultraviolet convergent: this 
has been achieved, as expected, by the substitution of R by Rd- However, many of the terms 
infrared divergent, and it is not at all apparent that such divergences will cancel between 
the various terms. To control these divergences, we evaluate the terms in dimensional 
regularization i.e. the p and k integrals are evaluated in a dimension d just above 2, and 
the resulting terms expanded in a series in (d — 2). We show below that while there are 
individual terms of order l/(d — 2) 2 and l/(d — 2), they do cancel among each other. 

By a series of elementary algebraic manipulations (including splitting apart some of the 
denominators by the method of partial fractions), all the terms in ( C15 ) can be expressed 
in terms of two basic integral expressions. These are 
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(C17) 



it is easy to see that both / and J are invariant under all permutations of their arguments 
i.e. I(a,b) = I(b,a), J(a,b,c) = J(c,a,b) = J(c,b,a) etc. In terms of / and J, the result 
( C15|) takes the form 
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dRD 3 dRD 

We can evaluate the needed values of I and J by standard methods which transform the 



momentum integrals into integrals over Feynman parameters Jl4 
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o 2/2 + 2/i(l — 2/i) (1 — 2/ 2 ) 



(C19) 



The last two expressions have been given directly in d = 2, as J(Rd, Rd, Rd) does not have 
any poles in (d — 2) and also does not appear in combinations multiplying poles in ( |(J18j ) . 
Inserting ( |C19| ) into ( |C18| ) and expanding in powers of (d — 2), we find that all double poles 
and poles in (d — 2) cancel, and we obtain our final result generalizing QUI Op : 



Ps{T) 
T 



3 

Gd~ 



Qd 
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(C20) 



3. n > 3 



The argument is more subtle for these cases with non-Abelian symmetry. We now ex- 
pect that at length scales ^> 1/V-R.d all longitudinal fluctuations will freeze out, and the 
transverse fluctuations will map onto those O(n) non-linear sigma model. This model has a 
dimensionless coupling constant g, and for g <C 1 there is a large correlation length of order 
(see e.g. Ref. 



iR 



D 



[n - 2 )g 
2tt 



1 V(n-2) 



2tt 



cxp 



(n-2) 9/ 



(C21) 



where we have chosen 1 / \/R~d as the natural short distance cutoff of the non-linear sigma 



model. For 1/a/Rd ^ \%\ ^ the two point- correlations behave as [13 

-, (n-l)/(n-2) 



($ Q (x)$ Q (0)) ex 



oc 



(ti-2)(71b( t 1 



1 - ^— ^ln( \x\JRd ) + ... 
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Let us compare this with the expression obtained from 



3T 

(*«(aO*a(0)) = ^ 



l-^^ln(NV^) + 



(C22) 

which yields in the same regime 

(C23) 



Comparing (|C22| ) and (|C23|) , we can obtain the missing prefactor in ( |C22|) , and also get 



Qd 
9 = T 



(C24) 
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Actually it is possible to do better, and actually fix the missing constant in ( U21|) pre 



cisely. To do this, as was shown by Hasenfratz and Niedermayer |33|], we need to carry 
out exactly the calculation of Section p 2| and obtain the H dependence of the free energy 
density for n > 3. Now the analog of the replacement ( |C11|) in the action (|1 . 7|) is 



(C25) 



a>2 



Next we will compute the H dependence of Z(H) in a perturbation theory in U, but will find 
that the structure of the answer is actually quite different from that found in ( |C15|) for n = 2. 
In the present situation one finds that all the infrared divergences in (1/H 2 ) lnZ(H)/Z(0) 
do not cancel, which correctly indicates that the renormalized stiffness p s (T) is strictly zero 
for all T > 0. Instead, the small H dependence of Z(H) is more complex, and we will show 
that \nZ(H)/Z(0) ~ H 2 \n(l/H). 

After substitution of ( |C25| ) into ( |1 . 7] ) , it is immediately apparent that at order U° Z(H) 
consists of two separate contributions. The first is the contribution of the $1,2 components 
and this is identical to that computed in Section |C~3 for n = 2, The second is the contribution 
of the remaining n — 2 components which, at this order, are simply free fields with 'mass' 
H 2 . So we get 



In 



2(H) 
Z(0) 



H 2 



3\R\ 



d 2 p 



TU J (2vr)V + 2| J R| 



[n 



+ (C26) 



(2tt 



where the first two terms are obtained by evaluating ( U14 ) to order U°, and the last term is 
the contribution of the remaining (n — 2) components. Now substituting Rd for R by using 
we obtain, as expected, an expression free of ultraviolet divergences: 



In 



Z(H) _ 3R D H 2 _ (w -2) r d 2 p 
Z(0) ~ 2TU 2 J (27y 



In 



'p 2 + H 2 

p 2 



H 2 



P 2 + Rd 



(C27) 



Finally, we evaluate the integral in the limit H 2 <C Rd, and express the result in terms of 
Qd using 



In 



Z{H) 
Z(0) 



H 2 



2Q 



3 (n - 2) 1 
+ — ^ln 



H 



D 



47T 



(eR D yi* 



+ ... 



(C28) 



We can now deduce the correlation length, £, from this result using the matching to the 
Bethe ansatz solution, as discussed in Ref |33|]: the result is 



1 fn-1 
' n - 2 



'R 



D 



e(n - 2)g 



1 l/(n-2) 



2tt 



exp 



(n-2)o 



(C29) 



where T(x) is the gamma function and g is given in 
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FIG. 1. Phase diagram of the theory Zq for d = 2, n = 1, 2 as a function of the temperature 
T and tuning parameter r. The quantum critical point is at T = 0, r = 0. The most important 
crossovers are represented by the dashed lines, and these occur at T ~ \r\ zv , where the dynamic 
exponent z = 1, and v is the correlation length exponent of the d + 1-dimensional classical ferro- 
magnet; these crossovers divide the phase diagrams into regions A, B, and C Region A is the high 
temperature of the continuum theory Zq, with the T —* oo limit taken after the short distance 
cutoff has been sent to zero to obtain the continuum limit; its properties are described by placing 
the r = scale-invariant critical theory at non-zero temperature. There are two low T regions, 
B, C, on either side of r = 0. The ground state for r > is a quantum paramagnet (or a Mott 
insulator, depending upon the physical system) with an energy gap; the dynamics in low T region 
C is described by a model of a dilute gas of thermally excited quasi- classical particles, and this 
shall not be discussed in this paper. The ground state for r < has long range order with {<p a ) ^ 
and the low T properties above it are described by a model of quasi- classical waves for n > 2 (for 
n = 1 a separate model of quasi-classical particles applies). There is a line of finite temperature 
phase transitions, T = T c (r), within region B at which the long-range order disappears; this is 
denoted by the full line. This T > transition is of the Kosterlitz-Thouless type for n = 2, and in 
the universality class of the two-dimensional classical Ising model for n = 1. The phase diagram for 
n > 3 differs only in that there is no line of T > phase transitions in region B i.e. T c (r) = 0, and 
long-range order is present only for T = 0, r < 0. The present paper uses the e = 3 — d expansion 
to develop a theory for the low frequency (u < T), long distance dynamics in region A directly in 
d = 2 for all n, using a model of quasi- classical waves (our model also contains the initial crossovers 
as T is lowered into regions B or C). In contrast, transport of the conserved 0(n > 2) charge in 
region A was discussed in II for small e; it was dominated by excitations with energy ~ T, and 
described by particles obeying a quantum transport equation. 
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FIG. 2. 'Tadpole' graph containing the only ultraviolet singularity of Z 
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FIG. 3. The 'dual' dimensionless coupling as a function of Q for n = 1 (full line), n = 
(dashed line), and n = 3 (dotted line). 
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0.5 1 1.5 

FIG. 4. Numerical results for (p s {T) /T) / (Qd /3) as a function of Gd/3. We used a square 
lattice of L x L sites with periodic boundary conditions and lattice spacing a. The lattice sizes 
used were L = 64 (stars), L = 128 (pluses) and L = 256 (X's). The dashed line is the locus of 
points where the Nelson-Kosterlitz jump ( |2.24|) is obeyed. The full line is the result of the small 
Qd expansion in ( |2.17| ). The filled circle indicates the position of the Kosterlitz-Thouless transition 
determined by the extrapolation to L = oo limit using the method described in the text and in 
Fig |||. In the approximation in which we assume that the T dependence of Qd is given by its 
leading value as T — > 0, Qd ^ 3T/p s (0), the scale on the horizontal axis becomes T/p s (0), while 
that on the vertical axis becomes p s (T)/ p s (0). 
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FIG. 5. Values of the coefficient A in ( |2.26| ) determined by fitting the L dependence of the 
measured p s (T)/T to ( |2.26| ). We fit the value of A to a linear function of Qd, and point where the 
line has the value A = 1 (indicated by the dotted lines), determines the position of the Kosterlitz 
Thouless transition. 
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FIG. 6. Scaling plot of the static susceptibility x(0) for n = 3 as function of Q. The dashed 
line indicates the prediction of weak-coupling expansion in (|2.1|) . 
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FIG. 9. As in Fig | but for n = 3 and Q = 20, 25, 30. For Q = 20, 30 we used two different values 
of Ra 2 (Ra 2 = 0.03, 0.04 for Q = 20, and Ra 2 = 0.04, 0.08 for Q = 30) and these are indicated by 
the presence of both dashed and full lines for these cases. The good agreement between the dashed 
and full lines is evidence that we are measuring the universal values in the continuum limit. 
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0.4 0.8 — 1.2 

CO 

FIG. 10. The Fourier transform of Fig |7| to frequencies, which yields the scaling function 
^ Sc(0,u,Q) in ( |1.26|) , where uJ = ujjc\fR. Results are for n = 1 and Q = 25,30,35,40. 
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